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^ . Using a parametrization of the Wilson loop with the minimal- area law, we calcu- 



o 

(N 



late the polarization operator of a valence gluon, which propagates in the confining 



^ ■ background. This enables us to obtain the infrared freezing (i.e. finiteness) of the 

running strong coupling in the confinement phase, as well as in the deconfinement 

(N : 

^ , phase up to the temperature of dimensional reduction. The momentum scale defining 

^-H ' the onset of freezing is found both analytically and numerically. The nonperturbative 

> : 

' contribution to the thrust variable, originating from the freezing, makes the value of 

CN ■ 

^ , this variable closer to the experimental one. 

o : 

■ 

o ■ 

■ I. INTRODUCTION 

^ ■ 

Qh! The path-integral representation for the Green function of a particle moving along a 

' closed trajectory and interacting with the gauge field yields the Wilson loop. When the 
> ■ 

•'-j . particle is confined by the gauge field, its Wilson loop obeys the area law. In this physically 
rS ' 

^ \ very important case the path integral cannot be calculated analytically, because finding 
the minimal surface for an arbitrary contour x^{t) m d> 2 dimensions is too complicated. 
Therefore, effective parametrizations of the minimal surface have been invented in the litera- 
ture (see e.g. , which means there exist certain physically motivated tricks to construct 
the minimal- area functional in terms of x^{t). 

To give an example we consider the following formula, which converts a double surface 
integral into line integrals with Stokes' theorem: 

/ da^u{x) / do^p{x)dydpD{x — x) = — f dx^ (p dXpD{x — x). (1) 
is Jt. Jc Jc 

Here D is an arbitrary function for which the integrals are finite, = and E is an 
arbitrary surface encircled by the contour C. As one can see, the choice D{x — x') = {x — x'Y 
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is the unique one, for which the derivatives on the L.H.S. of Eq. ([T|) are removed completely, 
and one obtains 

/ da^^{x) / da^y{x') = (h dx^ (t dx'^Xyx'^. 

Now, if C is a flat contour, then we can choose S flat as well, in which case the L.H.S. of this 
equation equals to 25'^, where S is the area of S. Therefore, for a flat contour, the minimal 
area reads 

S = \/ , where E^^ = / dx^^x^ (2) 



2 '^^ 

a" 



c 



is the so-called tensor area which is manifestly a functional of C only. Flat contours 
are a good approximation for particle trajectories when the particle is heavy. For example, 
the parametrization of Eq. ([2]) reproduces correctly the heavy-quark condensate js] and the 
mixed heavy-quark-gluon condensate ^\ in QCD. However, it certainly cannot be correct 
for a light particle, whose trajectory may deviate significantly from the flat one. 

A parametrization for light and even massless particles has been proposed in Ref. j^. 
Similarly to Eq. (I2l), it allows to express the area Smin of the minimal surface S in terms 
of a single integral. The main idea is to convert the proper time in the path integral to a 
length coordinate r G [0,-R]. After that, one can naturally parametrize S'min as an integral 
of the transverse direction, |r(r)|, along this coordinate: 

5'min = / rfr|r(r)|. (3) 

With such a parametrization of 5'min, we calculate in this paper the polarization operator 
of a gluon when it splits into two valence gluons. The latter move in a nonperturbative 
background, which confines them. In the absence of the confining background, the polariza- 
tion operator yields the standard Yang-Mills one-loop running coupling jsl. In the presence 
of the background, the running coupling goes to a constant in the infrared region, i.e. its 
logarithmic growth "freezes" [5|: 

47r An 

Here, 

3 

is the absolute value of the first coefficient of the Yang-Mills /3-function, m cx 
(string tension) is a nonperturbative mass parameter. 
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In this paper, we will show how freezing occurs and calculate the values of m and b both 
analytically and numerically. The paper is organized as follows. In section 2, we recall 
the derivation of ctsip"^) in the absence of the confining background, introducing the scalar 
polarization operator nfi.ee(p^)- In section 3, the world-line integrals for the polarization 
operator n(x, y) in the coordinate representation are evaluated. In section 4, freezing is 
discussed in detail. In section 5, we extend this approach to the analysis of freezing in 
the deconfined phase. In section 6, possible phenomenological consequences of freezing are 
discussed. Finally, the main results of the paper are summarized in Conclusions. 

II. POLYAKOV'S DERIVATION OF THE RUNNING STRONG COUPLING 

In this section, we recollect some steps of the derivation of based on the integration 

over quantum fluctuations of the Yang-Mills held (see also |^). The procedure starts with 
splitting the total Yang-Mills field into a background, A^, and a quantum fluctuation, 
a^, whose momentum is larger than that of the background. One can therefore substitute 
the Ansatz = + into the bare Yang-Mills action, 

So[A] = ^Jd'x{F^^AA]f, 

where -F'^i,[v4] is the QCD field-strength tensor and go the bare coupling. One separates 
into a slowly varying background field and fluctuations a°: 

A'', = A'; + a^. 

Fixing the so-called background Feynman gauge (D^a^)" = 0, one adds to the action the 
term 5g.f. = ^ J d^x [{D^a^yf, where (D^a^)" = d^al + f^'^Alal. The total gluon action 
reads 

5*0 + S'g.f, = 

= / d'^HKAA]? - 4a:{DMA]r - 2a; [S.^D^ + 2r'XM < + 0{a')} . 

To perform the one-loop renormalization of g"^, one has to integrate out the a^-gluons in the 
y4^-background. The full renormalized effective action can be written as 

s = So + s'- + s^- = J ^.j^f;MF;A-p). (5) 
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where the running coupling g{p) is the object of the calculation. The so-called diamagnetic 
part of the effective action for A^-fields has the form 



S'^''' = tr \n{-D^) = tr 
where we have used the decomposition 

-D^ = -d^ + A(i) + A(2) with A^^\x) = it" {d^Al + iA^d^) , /S.^'^\x) = {A^^t 



(6) 



aj.a\2 



(t") = —if . The paramagnetic part of the effective action reads 

S^-^ = ^tr In [1 + (-92)-^(2r^F^^)] = -^(-52);,^(2r^F^^(i/))(-92)-i(2/'^'^^F;',(x)). 

(7) 

The subscript "dia" describes the diamagnetic interaction of the A^-field with the orbital 
motion of the a^-gluons. This effect leads to the screening of charge and is present in 
the Abelian case as well. The subscript "para" is because this part of the effective action 
describes the paramagnetic interaction of the v4^-field with the spin of of the a^-gluons. It 
leads to the antiscreening of charge, which is a specific property of the non-Abelian gauge 
theories. The dia- and paramagnetic parts of the one-loop effective action can be written as 



s^^-^ = AT, . { ^^^^^ } . J ^F;.(p)F;.(-p)n,,ee 



ip') (8) 



where 



iifree(P ) - ^^^^ p2 " j (27r)4 q^{q + 



is the free scalar polarization operator. Equation ([5]) then yields 



I I h Kl 

— ---T7^.^-J^ (10) 



5(2 (p) 167r2 p 



where go = g{Ao). Introducing the renormalized cutoff A = Agexp ^— one finally 



arrives at the standard result as{p) 



III. CALCULATION OF THE POLARIZATION OPERATOR 



In reality, the gluon fluctuations a^'s do not appear as measurable excitations in the 
QCD spectrum and therefore must be self-confined, which means that the two a^-gluons 
propagating along the loop interact and form a colored bound state. This interaction has 
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Figure 1: Relative wavelengths of the three fields -B^, A'^, and a°. 



a one-gluon-exchange part plus a nonperturbative part, which may be related to string 
formation in the octet-octet channel coupled to a color octet. If one assumes Casimir scaling 
for the string tension, the corresponding string tension in the octet channel is related to 
that of the fundamental representation as (see e.g. [3|) cr = |crfund = 0.225 GeV^ with 
Cfund = 0.2 GeV^. The a^-gluons may be confined because of a stochastic background field 
B'^, whose momenta are even smaller than the momenta of the A^-gluons js] (see Fig. 1). 
The presence of this additional background field can be taken into account by substituting 
into the Yang-Mills action a modified Ansatz 



A; = B; + A'; + a';^. 

Accordingly, the definition of the one-loop effective action includes now the average over 
the background: S = — In^J Pa^e"'^''^])^, where (. . .)^ is some gauge- and 0(4)-invariant 
integration measure. In the course of the calculation of / Pa^e""^!"^', the presence of the 
background leads to the substitution d"^ D'^[B] in Eqs. ([6]) and (jTl), where (D^[i?]ai.)" = 
dfidl + f°'^'^B'^^a1. In the diagrammatic language, we are still considering the same one-loop 
diagrams, which contain only two external lines of the y4°-field, hut with the internal a^- 
loop receiving infinitely many contributions of the -B^-field. The latter appear through the 
path-integral representation of the operator (D^[i?])~^, which contains the Wilson line of 
the a°-gluon in the 5^-field. Such an a^-gluon is called valence gluon from now on. 
In the coordinate representation the scalar polarization operator has the general form 



nfi.ee(a;) 



(2^ 



Dq{x), where Dq{x) 



(11) 
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In an arbitrary background -B|J the polarization operator becomes 

The path-integral representation of this average reads (see e.g. jsl) 
{Il{x,y\B))^ = U{x,y) = J ds J ds J (VZf,)^y{Vz^)y^exp \- J dX— - J d\ 



JL 
4 



X 



X ( tr P exp 



I ( / d\z^Bl{z)e + / d\z^Bl{z)r 



(12) 



B 

Here, {Vz^)xy is the standard measure of integration over all paths z^{X) such that z{0) = y, 
z{s) = X. Namely, in d dimensions, 

{Vz,U = Jim n j J^^^ ^ zike), e = s/N. 

The S^-averaged term in Eq. (fT2l ) is the Wilson loop of the valence gluon. The confinement of 
the latter is refiected in the so-called area law. The Wilson loop (. . .)^ can be approximated 
by e"*^"^""", where 5*111111 = •S'mini^, z] is the area of the minimal surface encircled by the paths 
Zfj_{s) and ^^(s). Moreover, parametrization of S'min in the form of ([3]) leads finally to the 
restoration of the translation invariance of Il{x,y). 

To calculate Il{x,y) we introduce "center-of-mass" and relative coordinate of the glu- 
ons [sl: 

Ufj, = ^ _|_ _ ; ^fJ. = Zfj, — Zfj,. 

We define new integration variables, which have the dimension of mass: 

\x — y\ _ \x — y\ „. 
^=^^' ^ = ^^- (1^) 

Then the kinetic terms of the gluons become 
with 

R= \x — y\. 

Here, Ur = is the "reduced mass", and r is the distance to the point x along the 
line passing through x and y. Note that, unlike Schwinger's proper time s, which has the 
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dimension (length)^, the variable r has the dimension (length). The minimal area, 5'min, can 
be effectively parametrized with the coordinate r alone, after which the problem reduces to 
that of a Schrodinger equation with the potential cr|r|. The main new idea of the present 
paper is to use a path-integral approach for the calculation of n(x, y). To this end, we will 
apply the Cauchy-Schwarz inequality to Eq. ([3]): 



^ ■ 



1/2 



(14) 



Note that, in case of a (l+l)-dimensional classical-mechanics problem, this approximation 
works with a good accuracy. For two particles of mass m interacting through a linear 
potential, which move from {y,t) = (0,0) to the point (0, i?), one can write the equation 
of motion myi = —mi/2 = —o'. Substituting the solution, yi = —y2 = ^t{R — t) into the 
exact formula S'min = d^ivi — Z/2), we get S'min = Using instead our approximation. 



we obtain S„ 



< 



R 

R J dt{yi 




y2) 



1/2 



30m ' 



The relative error is therefore quite small, 



Encouraged by this observation, we return to the 4d case and eliminate the square root 
in the Cauchy-Schwarz inequality, Eq. (fT4l ). by introducing an integration over an auxiliary 
parameter A, which is sometimes called einbein. The expression for the polarization operator 
then becomes 



^ Jo H' Jo P 



dX 



TX 



exp 



-A 



/i + /i 



drii, 



R 



2 70 ^ 4A Jo 



(15) 



2 /- ^ 

Now, the integrals over u^{t) and r4(r) are free path integrals, while the integral over r(r) 
is that of a harmonic oscillator. We calculate these integrals by introducing the variables 



^ = aR?^"^ / ^/2fj,rX instead of A, a = fj,R/2 instead of fi, and b = fj,R/2 instead of fi. 



In terms of these variables 



U{x,y)^U{R) 



a 



167r9/2i?2 



f{aR' 



where 



f{<rR') 



d^ 



v/esinh^^/^^io 



dado \ I — — exp 
ab 



—a — b — 



aR'^\^ a + b 



2^ 



ab 



(16) 
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Figure 2: Integral f{aR ) from Eq. ([16 



When integrating analytically over a and 6, it has been found (cf. Appendix A) that, 
owing to the a ^ b symmetry, the corresponding saddle-point equations can be solved 
even when the pre-exponent is lifted to the exponent. However, the leading large- distance 
asymptotic behavior of n(i?) stems from the mere substitution of the saddle-point values 
without that lifting, a = 6 ~ into the pre-exponent (see Appendix A for details). This 
procedure yields 



a 



3/2 



-2<tRV? 



This integral can be evaluated by splitting the integration region into two parts: 



167r7/2i? 



(17) 



(18) 



where 



23/2 r^^-H_w 



Then, at cri?^ ^ 1, Ji = 0{e~~^'^^') is a subleading term. The integral I2 is saturated by its 
saddle point, ^ = 2R^ya/3, which lies inside the integration region. This yields 

.5/4 



m) 



a 



-2v^R ^i?2 ^ 1. 



(19) 



25/2 . 31/4 . ^3 . ^3/2 

Let us now consider the opposite limit of small distances. There, Eq. (A. 3) goes over to 
(2^)^^^' ^^'^ product of the three path integrals, Eqs. (A.1)-(A.3), yields 



/ 



{'Du^)xy(Vrf,)oo exp 



/i + 



7-2 r^r 

CLTU,, 

JO 



R ^2j^ r-R 

drrf, -— / drr'^ 



4A 
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/i/i 



(27ri?)2_ 



exp 



(/i + fi)R 



at aR^ < 1. 



Since this expression does not depend on A anymore, the remaining A-integration in Eq. (1151) 
results in a factor of 1. Equation (fTSl l then goes over to the free scalar polarization operator, 
Eq. (HI]): 



U{x,y) 



1 



(i/i / (i/i exp 



(/i + /i)i?' 



nfree(a; - v) at (T/?2 < 1. (20) 



Therefore, we have a formula for the polarization operator, interpolating between the lim- 
its (fT9D and (120D: 



n(i?) 



1 



-AvVR2 



with the following analytic values of the coefficients: A = 2y/3 ~ 3.46, B — '^^''^'^ 



(21) 



31/4 



6.75. 



IV. FREEZING OF THE RUNNING STRONG COUPLING AT ZERO 

TEMPERATURE 

The integral of Eq. (fT6l) has been also calculated with the Monte-Carlo integration routine 
Vegas jol in the interval 0.12 < aW' < 4.3. The interpolating curve is plotted in Fig. 2. The 
numerical fit to these data yields A = 3.38, which is very close to 3.46. The corresponding 
analytic and numerical values of the mass m in Eq. ([4]) are 

man = 2^3^ = 1.64 GeV, mr^nm = 3.38v^ = 1.60 GeV. (22) 

Note that the value of the analytically calculated freezing mass depends on the dimension- 
ality of space-time d as 

man = 2^/{d-l)a. (23) 

This is readily seen from the saddle-point of the integral I2 in Eq. (fT8l) by noticing that, in 
d dimensions, sinh^^^,^ — > sinh~ ^ in Eq. (flTl ). This result is a direct consequence of the 
Ansatz for the minimal area we use, Eq. f |T4l l. 

Freezing, as defined by Eq. dH), stems from the replacement of the free-gluon polarization 
operator, Eq. ([9]), by that of the valence gluon, 

nvai(p') = T7^1n^4^. (24) 
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Figure 3: The function F{q) (red online) and a fit to it (green online). 



Therefore, it looks instructive to compare the inverse Fourier image of this desired exact 
expression with our approximate result, Eq. (|2T1) . We have 

1"^ ds 

pA _|_ ^2 g 

ds .2 . 



d^pe-'P"" In ■ 



^^0 



d'^pe 



-ipx 



-e 4s 



K2{rn\x\), 



10 

where K2 is a Macdonald function. (When deriving the third formula in this chain, we have 
used the obvious fact that / d'^pe'^^^ = for x 7^ 0.) Therefore, Eq. (l24l) in the coordinate 
representation reads 



nvai(i?) 



m 



(25) 



The short-distance asymptotic limit of this formula coincides with Eq. (I20ll . As for the large- 
distance limit, we see that Eq. (l25ll has the same exponential fall-off as our result, Eq. (|2B . 
but a different pre-exponential i?-behavior. The ratio of Eq. (l2Tl) to Eq. (|25|1 at mR > 1 is 
oc ^/aR. However, at large distances in question, this discrepancy is unimportant. 

Finally, it is worthwhile to compare directly Eq. (i24l) with the Fourier image of Eq. (i2Ti l. 
Because of the logarithmic divergency, we compare the derivatives of the two expressions, 
which are UV-finite. Therefore, we compare '^'^^^^^ ^ = ~ le^ pi+m^ with where n(p^) 

is the Fourier image of n(_R): n(p^) = ^ dRR^Ji{\p\R)Il{R), J^'s are the Bessel func- 
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tions. Introducing the dimensionless variables x = \faR and q = \p\/y/a, one has 

dp^ a 

Fiq) = A-2 T-ll - ^2(^^)) - ^1 (1 + Bx'/') e-^^- (26) 

oTT^q^ Jq X l2 qx j 

This function has been calculated numerically for 0.02 GeV < |p| < 5 GeV, which corre- 
sponds to g G [0.042, 10.541]. Fitting the result by the function 



we obtain: di = 0.585, ^2 = 2.491. The interpolating curve for F{q), at the above-mentioned 
values of g, and the function (l27l) are plotted in Fig. 3. The value of m corresponding to the 
coefficient d2 is therefore 

m = rf2V^= 1.18 GeV. (28) 

It is closer to the phenomenological estimates (~ 1 GeV) {l^ than the values (p2l ). Further- 
more, the coefficient di defines a numerical prediction for the parameter b: 

b ~ 0.5856 = 6.435. (29) 

Therefore, the numerical analysis in the momentum representation yields an effective de- 
crease of b in the infrared region. For the final form of as{p) at zero temperature, we refer 
the reader to the end of section 5, where we discuss this result together with the one at 
T>T,. 

V. FREEZING OF THE RUNNING STRONG COUPLING IN THE GLUON 

PLASMA 



At temperatures above deconfinement, T > Tc, large spatial Wilson 



the area law. For pure gauge SU(3) Yang-Mills theory, T^. ~ 0.27 GeV [ll|. This behavior of 



loops still exhibit 



spatial Wilson loops is the main well established nonperturbative phenomenon at T > 



which is usually called "magnetic" or "spatial" confinement Of course, it does not 

contradict true deconfinement of a static quark-antiquark pair, since space-time Wilson 
loops indeed lose the exponential damping with the area for T > Tc. If points x and y are 
separated by a time-like interval, valence gluons in the polarization operator Il{x,y\T) are 
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not confined, and therefore as(p^) at > is given by the perturbative formula, without 
freezing. When x and y are separated by a space-like interval, magnetic confinement holds, 
and one expects freezing of as(p^) at < 0. 

We will start our analysis with the calculation of the spatial polarization operator of a 
valence gluon in the SU(3) pure Yang-Mills theory at temperatures higher than the tem- 
perature of dimensional reduction, T > Td.r. — 2Tc. QCD becomes a superrenormalizable 
theory in three spatial dimensions, where the renormalization of the dimensionful coupling, 
(73 = gVr, is exact in one loop. Apparently, this effective coupling is not asymptotically free 
and is unrelated to the study of freezing. The following calculation prepares the subsequent 
analysis of at Tc < T < Td.r. ■ 

Let us first consider the propagator of a free particle at temperature T from the origin 



to the point Rn 



in, Ri): 



(i?4 - ^nf 



R2 + (i?4 - 



i2n 



exp 



4s 



4s 

T I (X>z)Roexp 



s - 2 



—dX 

4 



(30) 



where ^ = Yl ■ Upon the Poisson resummation, one has 



1 



2T^/7^i 



exp 



(i?4 - Pn) 
4^ 



exp {—uj'^s + iuOnRi) , = 27mT. (31) 



When T — i> oo, only the zeroth term on the R.H.S. of Eq. fl3Tl) survives, which means 
dimensional reduction. The sum goes to 1, and 

T 



T 



ds I (Pz)Roexp ( - / —d\ 



s -2 



4nL' 



(32) 



where L = |R|. With the effect of magnetic confinement included, the polarization operator 
for T > Td.r. reads [cf. Eqs. (pi) and jllj)]: 



/•OO /•OO 

n(x, y|T) = T^ / ds ds I{s,s), where 
Jo Jo 



I(s,s] 



(•Dz)xy('Pz)yxexp 



4 



-dX 



-dX 



L 



drp- 



1/2 



(33) 



Here, as is the spatial string tension, whose ratio to the zero-temperature string tension, 
ao = 0.225 GeV^, is plotted in Fig. 4. The points x and y are spatially separated, y — x = R, 
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Figure 4: The curve interpolating the lattice data on the ratio of the spatial string tension to the 
zero-temperature one in the SU(3) quenched QCD as a function of T/Tc [inl- 



and p is a two-dimensional vector orthogonal to R. Two-dimensional vectors are denoted by 
an arrow, differently from three-dimensional vectors, which are boldfaced. This polarization 
operator can be calculated in a way similar to the zero-temperature one. Referring the 
reader for details to Appendix B, we present the final result: the polarization operator at 
T > Td.r. reads 

n(x,y|r)~^^e-'"^ ata,L2»l, (34) 



In Eq. (1341). 



n(x,y|T) 



T 
4^ 



at cXsL^ < 1. 



m 



(35) 



(36) 



i.e. Eq. (l23l) at li = 3 is reproduced. Equation (l35ll is nothing but the free scalar polarization 
operator, that is just the square of Eq. (I32ll . 

Let us now proceed to the physically more interesting range of temperatures, Tc < T < 
Td.r.. Here, like at T = 0, the exponential fall-off of n(x, y\T) (where x and y are four-vectors) 
due to magnetic confinement is relevant to the freezing of a^. The Euclidean path-integral 
representation for Il{x,y\T) at Tc < T < T^.r. can be constructed by using Eqs. (l30l) and 
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ds 



ds 



exp 



n,k 



(37) 



This quantity is calculated in Appendix B. At large distances, asL'^ S> 1, the result reads: 



Il{x,y\T) 



8V2tt^R^ 



(38) 



where R = y^Rl + L"^. Recalling that, in the physical Minkowski space-time, magnetic 
confinement holds only when i?^ is a space-like vector, one can always place the points x 
and y along some spatial axis, which makes L and R equal. The resulting formula for the 
polarization operator at large distances takes the form 



U{x,y\T) 



(39) 



Therefore, at the temperatures Tc < T < T^.r. and Minkowskian < 0, freezing of as(p^) 
takes place at the temperature-dependent momentum scale, which is analytically defined as 
m = l^los- The factor "2" under the square root in this formula is the number of spatial 
dimensions minus one, in accordance with Eq. ( l23l ). 

The limit of small distances is also discussed in Appendix B. At L = i?, the result reads 



n(x,i/|T) 



T 
4^ 



-\ 2 



coth{nTR) 



at a^R^ < 1. 



(40) 



(In particular, at T — oo, this result goes over to Eq. (I35l l. as it should do.) Since ^/oq = 
474 MeV, = 270 MeV, one can see from Fig. 4 that, at < T < Td.r., the condition 
(JsR^ <^ 1 automatically means also TR <^ 1. For this reason, at these temperatures. 



r. Therefore, the 



Eq. (I40l) can be approximated by its zero-temperature counterpart, p^^^^ip , 
formula for the polarization operator, which interpolates between this short-distance limit 
and the large-distance one, Eq. (I39ll . reads 



U{x,y\T) 



-2J2FIR 



(^1 + Tx\f2o,R^^ a.tTc<T < Td.r. 



(47r2i?2)2 

Note that this expression depends on temperature only implicitly, namely through o"s(T). 

The analysis of the polarization operator in the momentum representation can again 
be performed. Specifically, the factor (l + Bx^^'^) e~^^ in Eq. fl26l) should be replaced by 
(1 + TTy/2x^) e~^^^. Fitting the integral by the function (l27|l . we get the values di = 0.711, 
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Figure 5: The running coupling with freezing at < p < 3GeV for T = 0, T = Tc, T = 1.5Tc, and 
T = 2Tc. The curve at 3GeV < p < 4GeV is the experimental as{p), according to the web-site 



http://www-theory.lbl.gov/~ianh/alpha/alpha.html from Ref. 



d2 = 2.289. The freezing mass grows with temperature as the square root of crs(T): m(T) 



d2\/o's{T). When T varies from Tc to Td.r. — 2Tc, miT) varies from 1.09 GeV to 1.56 GeV. 
Finally, the parameter h ~ 0.7116 = 7.821 is larger than the one at zero-temperature, 
Eq. ([29]). 

Both at T = and Tc <T < Td.r., freezing modifies Eq. (fTOll as 



1 



m 



fl'^(p) fi'o 167r2 + 711? 
Defining the renormalized cutoff, A, through the bare one, Aq, as 

2n 



A = \ Aq + 7v? exp 



as(Ao)6 



one enforces as{p) 



Ait 



to take the value q;s(Ao) when p = Aq. It is then natural 



to choose a sufficiently large momentum scale Aq, where is practically unaffected by 
freezing and finite-temperature effects, and match it at that scale to the experimental value. 
Choosing Aq = 3 GeV, where as — 0.2524, we plot in Fig. 5 as{p) with freezing at T = 0, 
T = Tc, T = 1.5Tc, and T = 2Tc for p < 3 GeV, as well as the experimentally measured 



«s(p) [l7| for 3 GeV < p < 4 GeV. At at fixed p, one observes an increase of as{p) at T = Tc 
with respect to as{p) at T = 0, and a subsequent decrease with the growth of T. 
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VI. POSSIBLE PHENOMENOLOGICAL TESTS OF THE INFRARED 
FREEZING OF THE QCD COUPLING 

It has been pointed out in various papers [isl that the infrared behavior of the 
QCD running coupling can be used to estimate power-behaved corrections to various QCD 
observables. The essential parameters are few moments of the coupling in the infrared region. 
One class of the possible observables consists of event-shape variables like the thrust T: 

1 = max — - . 

Here the vectors p^'s are the momenta of the final-state hadrons and n is an arbitrary vector, 
which maximizes T. If the momenta of the hadrons form an almost collinear jet then, after 
the maximization, n will lie along the jet axis. The value of the above thrust variable becomes 
1 for an idealized pencil-like jet. Due to radiation of gluons, the observed T will be differing 
from 1. In perturbative QCD, these corrections from hard-gluon radiation at a certain scale 

n 

p can be calculated [16||. In addition to gluon radiation in the perturbative region, also 
soft-gluon radiation is present below a scale /iir ~ 3 GeV. Essentially the physics in this 
region is supposed to be parametrized by the freezing of a^, which we derived in this paper. 
Because of confinement, it seems to be impossible to map out this infrared running of the 
coupling exactly, i.e. for each momentum value. But measurements about the energy loss 
in fragmentation may allow to access an integral over the infrared region. The observable 
1 — T related to thrust T has an expansion in terms of the perturbative a^'^^^{p): 



1 - T 



= 0.334af + 1.02(af + O ((ar*(p))^) 

pert 



which acquires a hadronization correction due to soft gluon radiation 

1-T=1-T\ +—. (41) 

I pert p 

The nonperturbative higher-twist contribution A may be related to an integral over the 



infrared region of as{p) = z — "^2+^2 with freezing 



Jo 

where Cp = 4/3 is the Casimir operator of the fundamental representation of the group 
SU(3). Using for as{p) the parameters m and 6 at T = and T = Tc, we obtain 



A 



= 0.376 GeV, A 



= 0.395 GeV. (42) 

r=0.27GeV 
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In order to reproduce the experimental data on the energy losses through radiation in the 



0.5 GeV has been estimated 



T=0 



fragmentation process of c- and 6-quarks, a value A 

As a possible phenomenological application of freezing, let us evaluate the thrust of a 
two-jet event at the scale p = Mz. Using the value a^'^^^{Mz) = 0.12 one obtains the 
purely perturbative contribution 



T 



pert 



0.334af *(M^) + 1. 02 {a^/'\Mz)Y ~ 0.055, 



(43) 



which underestimates the experimental value measured at LEP 1 — T ~ 0.068. Inclu- 
sion of the higher- twist contribution according to Eq. f l4T! ). with the parameter A given by 
Eq. (1421) . yields for the full quantity 



T 



0.063, 



T=0 



(44) 



which is closer to the above-cited experimental value. 

If hadronization takes place in the quark-gluon plasma, then one can study the effect 
of a modified as{p) by investigating the energy loss due to soft radiation in the plasma. 



Accounting for the purely radiative higher-twist effect by means of A 
one gets at the M^-scale 



T=0.27 GeV 



from Eq. (142 



0.064. 



lT=0.27GeV 

If one instead uses as{p) from Refs. jl^, then one gets A 



(45) 



1.03 GeV, A 



T=0 



r=0.27GeV 



0.61 GeV, which overestimate the corresponding values (l42ll . as well as the experimental 



value A =0.5 GeV Hi 

T=0 

In conclusion of this section, although as{p) with freezing enhances the perturbative 
contribution ([43l) by 15%, the corresponding full values, Eqs. fl44l) and (l45l) . are still smaller 
than the experimental one. An increase of as{p) at T = T^, with respect to the zero- 
temperature case, enhances the energy loss in the quark-gluon plasma, which occurs due to 
the infrared radiation. However, the amount of this enhancement, quantified by (1 — T), is 
only 1.6%. 



VII. CONCLUSIONS 

This paper analyses the so-called IR freezing (i.e. finiteness) of the running strong cou- 
pling, in the confinement and deconfinement phases. The direct evaluation of the path 
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integral of a valence gluon confined by the stochastic background fields can be done by 
using parametrization ffT4l) for the minimal area. This parametrization reduces the path 
integral to that of the three-dimensional harmonic oscillator. In the deconfinement phase, 
at Tc < T < Td.r. — 2Tc, freezing of as{p) is present at Minkowskian < due to the 
so-called magnetic confinement. Upon the calculation of the path integral, we find the mo- 
mentum scales at which the freezing occurs in the confinement and deconfinement phases. 
Analytically, these scales are given by the same formula fl23ll . where d = 4 at T = and 
d = 3 at T > Tc (cf. Eq. (l36l) ). Numerically, the values of the freezing scales following from 
the fits in the momentum representation are smaller than the corresponding analytic ones 
and closer to the phenomenological value of 1 GeV (see Eq. ( l28l ) and the end of Section 5). 
The values of as(0) obtained at T = 0, Tc, l.STc, and 2Tc are 0.341, 0.377, 0.357, and 0.333, 
respectively. Full plots, which include the experimentally measured as{p), are presented in 
Fig. 5. Finally, we have estimated physical effect of the freezing on the thrust variable. The 
calculated nonperturbative contribution, which arises due to the soft radiation, brings the 
purely perturbative value of this quantity closer to the experimental one. 

Note in conclusion that one more application of the proposed method can be a calculation 
of mean sizes of the gluon bound states. Such bound states can be not only color-singlet 
(glueball), but also color-octet (quark-gluon) . The latter are suggested to be important 
ingredients of the quark-gluon plasma at temperatures <T < Tj , 



8|,l20fl. Work in this 



direction is now in progress. 
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Appendix A. Details of calculation of n(j;, y) at T = 0. 

The path integrals in Eq. (fTSll read as follows 



(/i + jj.)R 

r) ^^^L 

j (^^r4)ooexp (^-y^ drrl^ = 



(Al) 
{A.2) 
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y"(Pr)ooexp ^-y ^ drr 



4A 



R 



drr' 



UJ 



Att sinh ( ^ 



3/2 



(A.3) 



where u = a^y2fj,rR/ A is the frequency of the harmonic oscillator. Bringing these formulae 
together and passing from the integration over A to the integration over C, = ct-R^/^/a/^/vA, 
we arrive at Eq. (fTGl ) as it would look like before the introduction of the variables a and b: 



U{x,y) 



a 



211/2^9/2 v^sinh^/^^ 



djjidjji 



exp 



/i + /i 



R- 



a^R' 






2e 







{AA) 



It further turns out that the saddle-point integral over /i and n can be done even with 



the account for 1/v^Av in the pre-exponent. Indeed, we are dealing with the integral 

dfidfte''^^^''^^ , where 



fi + fi 



R 



a^R' 

2e 



lln 

2 



R \fi 



1 



The saddle-point equation §| = reads Rp,^ + — = 0. We can further use the 
fact that, since / is symmetric under n Ji, the saddle-point values of /x and Ji coincide. 
Therefore, setting in the last equation Jl = n, we arrive at a quadratic equation, whose 
solution is 

1 



^=AR 




{A.5) 



The sign " + " in front of the square root has been chosen because the integration region over 
fi and /2 is from to +oo. 

Further, as we are interested in the large-distance regime, aR^ > 1, we can consider 
separately the following three regions of integration over C, '■ 



^ < 1 < AaR^. At these values of ^, the saddle-point values (A. 5) are fi 



Accordingly, e" 



-/ 



/i ^ ^ 



— e Next, because / is a symmetric function of 

saddle— point ^ 

and /2, the pre-exponent of the saddle-point integral reads 



271 



2txo 



92/ 




dp? 


saddle— point 



3 i 



Altogether, 



dfid^e' 



2tt /2a 



8aR^ 



3 e 



-2(7ijV? 



(^.6) 
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Further, approximating ^^.^^3/2^ in Eq. (A. 4) by and changing the integration variable 



to t 



we have for the contribution to Il{x,y) from this region of ^'s 



n«(x,y) ~ — 



3/2 



where r{a, x) — dtt"" * is the incomplete Gamma-function. Using the known asymp- 
totics V{a,x) — > x"~^e~^ at x > 1, we finally obtain 



_l/2 -2f7R- / o 



327r7/2i?3 



{A.1) 



1 < ^ < AaE?. In this region of ^'s, Eq. (A. 6) still holds, but ^^^^^3^2^ in Eq. (A. 4) should 



be approximated by -^^^/2- This yields 



,^ — , where L = 
(327r7)V2i^ ' ^ 



at _M_22«i ^ o / 



3g 2<tH^ 
' 2 



(A8) 

Let us start with the analysis of Ji by representing the integration region as — 
Io° ~ lo ~ IlaR^- Then, the first of these three integrals can be done exactly and reads 



2Ko(2\'3aR), where here and below K^'s are the Macdonald functions. We can further take 
into account that the saddle point of e 2 « ^ which is ^ = 2y cr/Si?, lies between 1 and 
4(7 i?^. Owing to this fact, we can disregard y in the exponent of the integral , in the same 
way as we can disregard in the exponent of the integral J^pi2 ■ These approximations 
yield at aR"^ > 1 



h 2Ko{2\/3aR) - r(0, 2aR^) - r(0, 6aR^) 
Analogously, 

/9~ 



TT 



-2V3(tR 



3aR 



8aR^ 



ArJ-K,{2V3^R)- f d^e-^'^^'/^- r d^e 

V 3 Jo JactR'^ 



2aR^ 6aR^ 



-3^/2 



-2<7R2 



2e- 



-6(7i?2 



2aR^ 3 



at aR^ > 1. 



Equation (A. 8) then yields: 

nP)(x,y) 



a 



3/2 



(327r7)i/2i? 



1/2 



(3a)V4i?i/2 



2V3aR 
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5 e 



2aR^ y ' 8aR^ ) " 12 aR^ ' ^^'"^^ 

• 1 < 4(7 i?^ < ^. At these values of ^, the saddle-point values of // and /i, Eq. (A. 5), can be 
approximated as // = // ~ ^^j^. The exponent at the saddle point reads e~'^ 



. Again, owing to the symmetry of / under /i -(-^ /i, we find 



saddle— point 
saddle— point 



32^rR6- Therefore, the saddle-point integral over // and // reads 



-e 



Equation (A. 4) then yields for the contribution to the polarization operator, which comes 
about from this region of ^'s: 



VTT'e J4< 



4.i^2 ^V2^ 



The value of the saddle point of the exponent, ^ — -^^rji — , is smaller than AaR , for which 
reason we approximately have 



4 roo 



— "2" 



277^6^ ^ V 2' y 48V7r7ei?3 



(.1.10) 

Bringing together Eqs. (A. 7), (A. 9), and (A. 10), we can write the final result for the polar- 
ization operator n(a;, y) — H^^^ (x, y) + II^^^ (x, y) + Yi^^^ (x, y) : 



Il{x,y) 
where 



a 



5/4 



25/2.31/4.7^3.^3/2 I ^ + 4i^^' 



1 _ 23/2 ^1/2 

327r7/2 i?3 



e-V2 _ 5 . 2-1/2 ^1/2 



487rV2 



i?3 



Note that the correction O ^'^^^^ is negative. The leading term, O {e 
about from Eq. (A. 9). 



comes 



Appendix B. Details of calculation of n(x,?/) at T > Tp. 



Let us start with the polarization operator at T > Td.r.. Choosing for concreteness 
R = (L, 0,0), we define the relative coordinate r(T) = z(t) — z(t) = (ri(T),p(T)) and the 
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"center-of-mass" coordinate u = ^f^. After changing the variables (s, s) to yu) according 
to Eq. (fT3l) . we have for polarization operator (l33l) [cf. Eq. (fTSl) ]: 

U{^,y\T)^i^^y J^^^J^^pi^^,^), where (5.1) 



J(/i,/i) = / (r'u)xy(r'r)oo / ^^exp 



-A- 



vttA 

Carrying out the path integrations over r and u, we have 



+ / 2 



2 Jo 4A Jq 



drp' 



dX 



(5.2) 



Changing the integration variable A ^ ^ = Usi/ one can rewrite this equation as 



4V2^ 



--PP^/Jl + JjL 



^sinh^ 



exp 



2" 



/i + /i + 



1 1 
- + - 



(5.3) 



One should now substitute this expression into Eq. (B.l) and perform approximate saddle- 
point integrations over /i and Ji. This procedure means the insertion of the saddle-point 
values into the pre-exponent. In this way, one arrives at the following counterpart of Eq. (fTTI) 
at T > Td.r.: 



n(x,y|T) 



(5.4) 



871^ /^L Jo V^sinh^ 

d-l 

Note that, in three dimensions, sinh 2 ^ = sinh^ as explained after Eq. ( 1231 1. By splitting 
the integration region as in Eq. ( ITSl l. the last integral can be approximated as I1 + 12, where 



d^ 



To evaluate these integrals notice that, at temperatures T > T^ r,, which we are currently 
considering, the spatial string tension is parametrically 0"^ oc (73 ~ T^. Numerically, the ratio 
as/(Jo is larger than 2 at T > Td.r. (cf. Fig. 4). Next, for magnetic confinement to hold, 
the spatial Wilson loop should be sufficiently large, in particular the distance L should be 
> 1 fm. For such distances, 2a sL'^ > 'iaoL? ^ 1, and we have 
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Bringing these expressions together and neglecting the subleading terms ~ ^ (mLp'' i 



arrive at Eq. (I34l l. 

In the opposite limit of small distances, Eq. (B.2) yields 

Inserting this expression into Eq. (B.l), we have 



(5.5) 



n(x,y|T) 



^2 



327r3L 



(i/i /"°° (i/i 



(/i + /x)L 



at CToL^ < 1. 



Straightforward integrations over jj, and /i in this formula lead to Eq. (I35ll . 

Let us now consider the temperature interval < T < T^ r.. Changing in Eq. (j37l) the 
variables (s, s) to (/i, /i), we have 



n(x,y|T) = — 



n.k 



Ji) exp \ -— [/i(i?4 - f3nf + /i(i?4 - [3kf] [ . 



(S.6) 



Using for /(/i, /i) representation (B.3) and performing the saddle-point integrations over /x, 
/2 as above, we arrive at the following analogue of Eq. (B.4): 



n(x,i/|T) 



3/2 
(Ts 



16n^V2^LJo ^3/=^sinh^ 



n,k 



e « 



where 



1 + 



-R4 — /?n 



When splitting the ^-integral as in Eq. (fTSll . 



^3/2 gin]^ ^ 



e « 



^5/2 



e « 



dj 
1 



-((/3„+(/3fc) 



(5.7) 



we notice that, at temperatures Tc < T < Td.r. of interest, the ratio cTs/cto ranges between 1 
and 2 (see Fig. 4). Therefore, 

asL'^i^n + Vk) > 2a sL^ > 2aQL^ > 1 at L > 1 fm. 

For this reason, the second of the two integrals on the R.H.S. of Eq. (B.7) is again saturated 
by its saddle point, ^s.p. = ^JasL'^i'^n + V^fc), and yields the leading exponential fall-off with 
L, whereas the first integral yields a subleading Gaussian term. Disregarding the latter, we 
have for the leading exponential fall-off: 

^-2^ asL^{iPn+Vk) 



n(x,i/|T) 



E 
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Now, as we have just seen, the argument of the exponent here is much larger than 1. 
Therefore, one can safely restrict oneself to the y^o-terms in the sums over n and k. Since 
V^o = R/L, where R = \/R^i + L'^, we arrive at Eq. ( l38l ). 

In the opposite case of small distances, inserting Eq. (B.5) into Eq. (B.6), we have 



U{x,y\T) 



1 



647r4L2 



n.k 



-T 



1 



d/xexp 



167r4 ^ [L2 + (i?4 - (3nY] [L^ + (i?4 - (3ky 



at cr,i?^ < 1. 



(5.8) 



Therefore, we have recovered in the short-distance limit the product of two free scalar 
thermal propagators. Note that the sums in Eq. (B.8) can be done analytically. Setting 
L = R, Ri = 0, we obtain Eq. (1401) . 
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